AD-A204  300 


ffnn  Ffi  E  copy 


REPORT  SRL-P-01-19S9 

COMPACT,  LIGHTWEIGHT  CO2  LASERS 
FOR  SDIO  APPLICATIONS 


Prepared  by 
Dr.  Jonah  Jacob 


SCIENCE  RESEARCH  LABORATORY,  INC. 
15  Ward  Street 
Somerville,  MA  02143 


31  January  1989 


PROGRESS  TECHNICAL  REPORT 
Period  for  September  30,  1988  to  January  31,  1989 
Contract  Number  N00014-88-C-0727 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


Research  Supported  by 
STRATEGIC  DEFENSE  INNIATIVE/ 
SCIENCE  k  TECHNOLOGY 


Managed  by 

OFFICE  OF  NAVAL  RESEARCH 


DTIC 


ELECTEI 
FEB  1  7198^ 


H  ^ 


“The  views  and  conclusions  contained  in  this  document  are  those  of  the  authors  and  should  not  be  inter¬ 
preted  as  representing  the  official  policies,  either  e.xpressed  or  implied,  of  the  Strategic  Defense  Initiative 
Organization  or  the  U.  S.  Government," 


89  2  16  188 


RF  STABILIZATION  OF  MOLECULAR  DISCHARGES 


"^wo  rate  equations  axe  important  in  determining  the  stability  of  the  discharge.  The 
first  is  the  electron  production  and  loss 

2/3  /I  \ 

=  urimTie  -  an^  -  PUe  (1) 

and  the  second  is  the  metastable  production  and  loss 

=  {(Tv)ntna  -  n,n/rm  (2) 

where  rig  is  the  electron  density,  n„,  is  the  metastable  density, is  the  recombination  rate, 

0  is  the_attachment  rate,  \crvj  is  the  electron  impact  metastable  production  rates  constant 

'H-  Cl  u*-  r. 

and  T„j  is  the  metastable  lifetime.  The  ionization  rate  constant  v  in  Eq.  (1)  is  assumed 
to  be  the  result  of  metastable  ionization.  Since  the  metastable  levels  have  a  much  smaller 
ionization  energy  than  the  ground  state  the  ionization  rate  of  the  discharge  is  dominated 
by  electron  impact  ionization  of  the  metastables.  The  stability  of  Eqs.  (1)  and  (2)  will  be 
analyzed  assuming  that  the  applied  discharge  is  sinusoidally  varying.  Hence  the  various 
rates  are  determined  by  the  local  RF  electric  field  of  frequency  w  which  is  generated 
externally  by  the  system  shown  in  Fig.  1.  From  the  analysis  performed  in  Phase  I  Science 
Research  Laboratory  (SRL)  has  shown  that  a  current  source  is  volumetrically  stable  and 
so  the  initial  analysis  will  assume  an  RF  current  source,  f'  •  '  v  c  y  4  <■  '  '  ^ 

We  shall  assume  now  that  the  rates  describing  the  variation  of  rig  and  in  Eqs.  (1) 
and  (2)  (i.e.  arig,  3,  -frif/nm)  are  much  smaller  than  u;/27r,  so  that  one  can  use  time 

averaged  rate  constants  in  these  equations.  On  the  other  hand,  we  shall  assume  that 
u}/2it  is  much  smaller  than  the  electron  energy  distribution  equilibration  rate  ii'inciasUc)i 
and  therefore  in  the  above  mentioned  time  averaging  process  one  can  use  the  rates  corre¬ 
sponding  to  the  stationary  discharge  with  the  electric  field  equal  everywhere  to  the  field 
in  the  RF  discharge  at  a  given  time.  Thus  we  assume 

<<  ^'me/aatic  (3) 
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A  ^  ; 
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Assuming  a  100  Torr  discharge,  urim  ~  10®  sec"*  and  =  10**  sec"*,  so  that  1  MHz 
«  f  =  u;/2Tr  «  102  GHz. 

At  this  stage  it  will  also  be  assumed  that  the  vacuum  wavelength  of  our  oscillator 
(A  =  27rc/u;)  is  much  larger  than  L,  so  that  in  the  vicinity  of  the  discharge  one  can 
neglect  the  relativistic  retardation  of  time  due  to  the  extended  voltune,  and  treat  the 
electromagnetic  problem  as  being  quasistationary.  For  Z  =  1  m,  the  condition  L  «  X 
yields  /  <<  c/L  =  3  •  10*  sec"*.  So  that  we  are  limiting  our  discussion  to  frequencies  in 
the  range  10  MHz  <  /  <  100  MHz. 

The  quasistationarity  condition  means  that  in  the  Maxwell  equation  for  V  x  B  one  can 
neglect  the  displacement  current  and  use  the  "‘magnetstatic”  law 

=  (4) 

while  E  is  given  by  the  Faraday  law  and  the  quasi-neutrality 


Vx  £■= 


IdB 

c  dt 


V-  E=0 


(5) 

(6) 


Equations  (4),  (5)  and  (6)  yield 


72-^  j 


E= 


c*  dt 


(7) 


where  the  current  density  in  the  plasma  is  given  by  j  =  a  E,  with  the  conductivity  cr  = 
—eritfi.  The  simplest  model  for  the  mobility  ^  will  be  used  to  give 


a  = 


e^Tig 

- ^  S) 

mvcou 

Now  we  shall  consider  the  steady  state  of  Eqs.  ( 1-2)  ass\iming  cylindrical  symmetry 
and  that 

E=  Re{Eot-^ye, 


Therefore,  all  the  parameters  of  the  steady  state  (7o,Bo)  are  functions  of  r  alone. 

In  steady  state  one  obtains 

=  irn^o  (9) 
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— 


U^T  —  Ct 

All  the  rate  constants  are  temporally  averaged  over  the  voltage  oscillation. 
Returning  back  to  Eq.  (7),  the  steady  state  equation  gives 


(10) 


\  d  f  dE. 


dEo\  .47re^n* 

dr  )  mc^ 


UJ 


r  dr  \  dr  /  me*  u^oii 
The  solution  of  this  equation  is  (assuming  ~  const(r)) 


where 


Eg  =  AJoixr) 


2  _  .47re^neo  lj  _  2i 


X  =  I 


and  8  is  the  collisicnal  skin-depth 


VcoU  8^ 


{  cmvcoii  \ 


^  ^  ^  /  "''■"■'CO**  \ 

V^27rcra;  \2ire^neoi^ ) 


(11) 


(12) 


(13) 


(14) 


The  constant  A  in  Eq.  (12)  can  be  foimd  from  the  condition  that  the  magnetic  field  on 
the  surface  of  the  pleisma  equals  to  2/o/ci?,  where  Ig  =  Vo/Z  is  the  current  given  by  the 
external  current  source.  From  the  Faraday  law 


so  that 


and  therfore 


zc 

= - xAJiixr) 


A  =  i 


2IoU> 


c'^BxJiixB) 


In  the  limit  of  low  frequencies,  when 


i?  c^mucoii 

8  2neWoB? 


(15) 


(16) 


(17) 


(18) 


we  can  use  the  small  argument  expansion  of  the  Bessel  function  to  obtain  the  expression 
for  Eg: 


[-1(9 


(19) 
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Condition  (18)  is  important,  since,  as  follows  from  (19),  it  allows  one  to  have  a  radially 
uniform  discharge  in  the  steady  state.  For  example,  in  the  case  of  a  100  Torr  laser  discharge 
{I'coli  ~  3-  10^^  sec~^),  R  =  Z  cm  and  rieg  =  10‘^  cm~^,  (18)  yields  w  <<  2-10^®  (/  <<  3 
GHz)  which  is  within  the  quasistationary  limit  (/  <  100  MHz)  described  ezurlier. 

When  (18)  is  satisfied,  Eq.  (19)  yields 


Eo(r)  ~ 


\RJ 


const (r) 


It  is  precisely  this  inverse  dependence  of  Eg  on  n*  which  assures  the  global  stability  of 
the  discharge  in  the  current  source  scheme.  Any  increase  (or  decrease)  of  the  total  density 
leads  to  the  decrease  (or  increase)  of  Eg  and  therfore  to  the  decrease  (or  increase)  of  the 
electron  and  metastable  production  rates.  This  suppresses  the  initial  perturbation.  The 
details  of  such  a  global  stabilization  by  the  current  source  are  similar  to  those  already 
considered  previously  for  the  DC  discharge.  In  contrast  to  the  DC  discharge,  however, 
the  RF  current  source  also  allows  to  control  the  local,  streauner  type  instabilities  via  the 
skin-effect.  This  analysis  will  be  presented  next. 

1.1  Stability  Analysis 

Assume  that  the  lairge  skin-depth  condition  (18)  is  satisfied  for  a  given  steady  state 
so  that  Hgo,  rijnoi  O’oi  Eo  “almost”  (to  order  {R/S)'^)  uniform  throughout  the  plasma 
volume.  One  can  introduce  a  small  local  streamer  type  density  perturbation  of  a  char¬ 
acteristic  transverse  size  A.  The  diffusion  term  will  be  added  to  the  rate  equations  since 
non- uniformities  axe  of  interest.  Hence  one  can  write 


=  irimrit.  -  anl  -  +  Da'^^n^ 


“1" 

at  T 


where  Dg  and  Dm  are  the  ambipolar  and  met2istable  atomic  diffusion  coefficients  respec¬ 
tively.  Dm  is  of  the  order  ~  10^  cm^/sec/Torr  while  Dg  may  probably  be  10  times 

larger  {Dg  ~  10®  cm®/sec  Torr). 
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The  perturabation  in  n*  leads  to  a  local  conductivity  perturbation  which  in  turn  results 
in  the  decrease  or  increase  in  the  (rms)  electric  field  value,  depending  on  the  sign  of  the 
density  perturbation.  The  effect  is  similar  to  the  skin-effect  and  can  be  used  to  suppress 
a  possible  absolute  instability  in  the  discharge  (as  that  associated  with  the  attachment). 
This  stabilization  phenomenon  will  be  considered  in  more  details.  In  what  follows,  the 
diffusion  of  metastables  will  be  neglected  for  simplicity.  It  will  also  be  assumed  that 
is  the  fastest  rate  in  Eqs.  (20  and  (21),  i.e., 

i  »  UTlrrxo,  ""eo,  0,  I?a/A^  (22) 

Typically  the  rate  is  in  the  few  MHz  region,  while  the  other  rates  belong  to  the 
sub-MHz  scale.  For  times  t  >  l/zs,  condition  (22)  Eq.  (21)  to  be  replaced  by 


nm{t)  ^ 


(23) 


(this  is  the  fast  scale  analysis  of  the  previous  study)  Then  Eq.  (20)  becomes 

^  =  (t-7r  -  a)nl  -  /Jn*  +  Da^'^n^ 


(24) 


If  rit  =  const,  Eq.  (24)  again  yields  the  steady  state  Eq.  (10). 

At  this  point  the  perturbation  analysis,  will  be  employed  i.e.  solution  of  Eqs.  (7)  and 
(24)  in  the  form 


Tlf  —  rieo  4"  1 

i=(E.  +  E,)e-‘eJ 

where  all  the  perturbations  (ndjEi)  vary  slowly  in  time 


<<  1,^  = 


din  ^ 

\u;l 

dt 

The  perturbed  equations  become  (after  the  linearization): 


=  -i 


47re^u; 

c^mvcoii 


(n^lEo  +  ritoEi)  + 


4n-e^  /  dn^i  dEi 

c^mucou  V  °  dt 


Origi 

&t 


(25) 


(26) 


=  [2(i/7r  -  Q)n„  -  /?]  -f  [{v-)t  -  a)'n\^  -  0'n^o]  Ei  (27) 
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Here  (...)^  =  dfdEo  (...)•  Finally,  neglecting  the  selected  term  in  Eq.  (26)  and  using  the 
steady  state  relation  (10)  in  rewriting  Eq.  (27)  one  obtains 

+  (29) 

The  skin-depth  S  hzis  again  been  introduced  in  Eq.  (28). 

At  this  stage  a  solution  of  Eqs.  (28)-(29)  in  the  form  ne],jEi  ~  e**  ’’  axe  searched  for, 
where  \k\^  is  of  0{‘lt^jh.)  Then  one  has 


,  2  .  Cl  •“*  -^o 

^  "b  El  —  Ugi 


where 


^  =  (/3  -  D,P)n,i  4-  l3n,oB'Ei 


B'  =  -  a)'  _  ^ 


{u‘yT  —  a)  /? 

By  substituting  Ei  from  (30)  into  (31)  one  obtains 


7-17  2  4:/3neoEoB\  ,  .2l3nt„EoB'k^6'^ , 

■rfT  =  ‘  -  w'+4  )  + 

Since  k  >  2-k/R  and  S  »  R  one  finds  (ka)*  >>  4  and  therefore  (33)  reduces  to 


_  n  ,2  ^dn^oEpB'  2l3n^oE„B' 

dt  ^  ^  ^ 

The  last  equation  shows  that  a  stable  operation  corresponds  to  the  case 


,  2  4l3nfoE„B' 


^  H  ^  -  D.k^  -  <  0  (35) 

One  observes  that  there  exist  two  stabilizing  phenomena  in  the  system.  The  first  is 
due  to  the  diffusion  and  is  effective  in  stabilizing  short  scale  density  perturbations  The 
diffusion  simply  smoothes  the  short  scale  non-uniformities.  The  second  phenomenum  is 
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basically  the  skin  effect.  It  is  more  important  in  suppressing  larger  scale  perturbations, 
espesizdly  those  of  size  comparable  to  the  radious  of  the  plasma  column. 

The  interplay  between  the  two  stabilizing  processes  can  be  demonstrated  by  drawing 
the  curve  showing  the  dependence  of  $  =  (see  Eq.  (35)).  A  typical  situation  is 

shown  in  Fig.  2. 

$(fc^)  reaches  its  maximum  value 


at 


,2  _  (S0n,oEoB'V>^ 

V  ) 


1/3 


(36) 


kL  = 


(37) 


Therefore,  in  case  km  >  2'it[R^  the  stability  condition  for  the  entire  k  -  space  becomes 


(38) 


Note  that  the  stable  operation  is  easier  at  higher  frequencies  and  pleisma  densities. 
Note  also  that  the  final  value  of  Da  is  important,  since  (38)  can  not  be  satisified  for 


0. 
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